Evidence of fractional matching states in nanoperforated Nb thin 
film grown on porous silicon 
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PACS 74 . 25 . Fy - Transport properties of superconductors 
PACS 74.25.Qt - Vortex lattices, fiux pinning, flux creep 
PACS 81 . 05 . Rm - Porous materials 

Abstract. - Resistive transitions have been measured on a perforated Nb thin film with a lattice 
of holes with the period of the order of ten nanometers. Bumps in the dR/dH versus H curves 
have been observed at the first matching field and its fractional values, 1/4, 1/9 and 1/16. This 
effect has been related to different vortex lattice configurations made available by the underlying 
lattice of holes. 



Introduction. — The nuclcation of Abrikosov vor- 
tices [1] in the mixed state of type-II superconductors with 
periodic artificial pinning centers attracted a great atten- 
tion since 1970's. Recent progress in the fabrication of 
nanostructures provides the possibiUty to reahze super- 
conducting thin films containing artificial defects as pin- 
ning sites with well-defined size, geometry and spatial ar- 
rangement [2,3]. Vortex pinning was extensively explored 
by many groups to develop a fundamental understanding 
of flux dynamics and for its relevance in applications which 
require enhancements of the critical current density. Thus, 
several types of artificial pinning centers, such as square, 
rectangular or triangular arrays, have been introduced in 
a controlled way in the superconducting films. In particu- 
lar, the use of regular array of pinning centers such as an- 
tidots [2,4-6] or magnetic dots [3,7-9] brings to new com- 
mensurability effects, which give additional insight into 
the pinning properties of vortices. The most notable phe- 
nomenon for these studies is the so-called matching ef- 
fect which occurs when the vortex lattice is commensurate 
with the periodic pinning array. This situation occurs, in 
particular, at fractional or integer values of the so-called 
first matching field Hi = $o/ao' when the applied 
field H corresponds to one flux quantum, $0 = h/2e, per 
unit cell area, Oq, of the pinning array. Here ap is the lat- 
tice constant of the pinning arrangement. As a result, at 



the matching field, the critical current density, Jc, is dras- 
tically enhanced [3,4,10] and moreover, as a consequence 
of the Little-Parks effect [11], the upper critical magnetic 
field is increased at the matching values. Recently anti- 
dot arrangements with a big variety of symmetries have 
been investigated. Matching effects have been reported in 
perforated Nb thin films for antidots lattices with short 
range order [12], or quasiperiodic fivefold Penrose struc- 
tures [13]. Moreover asymmetric pinning arrays have been 
suggested as superconducting rectifiers [14]. 
If the artificial structure of defects is created by litho- 
graphic technique, the matching fields are usually in the 
range of a few oersteds. For this reason, matching ef- 
fects are observed in a very narrow temperature region, 
close to the critical temperature Tc, for a reduced value 
t — T/Tc > 0.95. In order to both increase the match- 
ing field and decrease the temperature where the effect 
is present, the period of the pinning structure should be 
reduced to less than 100 nm. This gives, in fact, the possi- 
bility to increase Hi up to 1 tesla or even higher. A reason- 
able method to achieve this goal is to use self-assembled 
substrates, such as, for example, AI2O3 templates with 
characteristic features in the nanometric scale [15]. The 
pore diameter in AI2O3 substrates could easily be varied 
in the range 25-200 nm with porosity (i.e. interpore spac- 
ing) around 50%, and this gives the possibility to achieve 
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matching fields of tliousands of oersteds [15]. To prepare 
AI2O3 substrates bulk Al [16], Al foils [17], and deposited 
thick Al films were used [15]. 

Very recently, another very promising material for self- 
assembled substrates and an optimum candidate for the 
Nb growth was proposed, namely, porous silicon (PS) [18]. 
PS is constituted by a network of pores immersed in a 
nano crystalline matrix [19] and it is a material which of- 
fers a considerable technological interest in different fields, 
as for instance micro and optoelectronics [20] and gas sens- 
ing [21,22]. The diameter of pores, 0, in PS can easily 
be varied from 200 nm down to 5 nm by using substrates 
with appropriate doping (n or p) and different regimes of 
anodization. The porosity, in fact, can be varied in the 
range 30-90% by adjusting parameters such as the acid 
solution, the anodizing current density and the illumina- 
tion of the substrate during the anodization. The regu- 
larity of the pores arrangement, however, is of the order 
of 10% lower than the one observed in AI2O3 templates 
obtained by electrochemical oxidation [23]. It has been 
demonstrated [18] that thin Nb films deposited on PS sub- 
strates can inherit their structure. The resulting samples 
then consist of porous Nb thin films with in plane geo- 
metrical dimensions, oq and 0, comparable with the su- 
perconducting coherence length, £,{T). In these samples, 
matching fields of the order of 1 Tesla were experimentally 
observed [18]. 

Aim of this work is to deepen the study of the matching 
efltect in superconducting Nb thin films deposited on PS. 
Superconducting properties were investigated by trans- 
port measurements in the presence of magnetic fields ap- 
plied perpendicularly to the samples surface, down to t = 
0.52. As a consequence of the high density of the pore 
network, the (H,T) phase diagram presents a deviation 
from the classic linear dependence. This effect appears 
at the matching field Hi ^ 1 Tesla, a value larger than 
those typical of periodic pinning arrays obtained both by 
lithographic techniques and by using another kind of self- 
organized templates. Moreover a new effect related to the 
commensurability between the vortex lattice and the un- 
derlaying pinning structure was found. It consists in the 
appearance of pronounced structures in the derivative of 
the R{H) curves, dR/dH, which can be observed in cor- 
respondence of the first matching field and its fractional 
values. 

Fabrication. — Porous layers were fabricated by 
electrochemical anodic etching of n-type, 0.01 flcm, 
monocrystalline silicon wafers. The electrochemical dis- 
solution was performed in 48% water solution of HF, ap- 
plying a current density of 20 mA/cm^. The anodization 
time was chosen in the range of 0.5 - 4 min in order to 
get porous layers with a thickness ranging from 0.5 to 
4 /xm. The pores extend on a surface of about 1 cm^. 
The integral porosity was estimated by gravimetry to be 
of about 50% [24]. The resulting porous substrates have 
0=10 nm and gq = 40 nm. For this lattice, if the formula 



Hi — ^o/ttQ for the square lattice is used, the expected 
first matching field is Hi — 1.3 Tesla. 
Nb thin films were grown on top of the porous Si sub- 
strates in a UHV dc diode magnetron sputtering system 
with a base pressure in the low 10^^ mbar regime and 
sputtering Argon pressure of 3.5 x lO"'^ mbar. In order 
to reduce the possible contamination of the porous tem- 
plates, the substrates were heated at 120°C for one hour 
in the UHV chamber. The deposition was then realized 
at room temperature after the cool off of the substrates. 
Films were deposited at typical rates of 0.33 nm/s, con- 
trolled by a quartz crystal monitor calibrated by low-angle 
reflectivity measurements. Since the effect of the periodic 
template would be reduced when the film thickness, d^h, 
exceeds the pore diameter, 0, [18] the Nb thickness was 
chosen to be 8.5 nm for the sample analyzed in this paper. 
A reference Nb thin film of the same thickness was grown 
on a non-porous Si substrate in the same deposition run. 

Experimental results and discussion. The su- 
perconducting properties were resistively measured in a 
^He cryostat using a standard dc four-probe technique on 
unstructured samples. The critical temperature was de- 
fined at the midpoint of the R(T) transition curves. The 
value of the transition temperatures of the film grown on 
the porous substrate and of the reference sample in the ab- 
sence of the magnetic field were Tc = 3.83 K and Tc = 4.53 
K, respectively. The critical temperature depression in the 
case of the porous sample is consistent with what already 
reported in literature for films grown both on AI2O3 [15] 
and on PS [18]. The first step for the characterization 
of the behavior of the porous Nb sample in the presence 
of perpendicular magnetic field is the determination of its 
{H,T) phase diagram. The temperature dependence of 
the perpendicular upper critical field, Hc2±, was obtained 
performing resistance vs. field, R{H), measurements at 
fixed values of the temperature with a temperature stabil- 
ity of 1 mK. Hc2i_ was defined at the midpoint of each of 
the R{H) curves. 

In Fig. [T]the {H,T) phase diagrams of the Nb thin films 
are shown. In general, the perpendicular upper critical 
field of superconducting films of thickness d obeys a lin- 
ear temperature dependence, Hc2±{T) ~ (<I'o/27r^Q|j)(l- 
T /Tc) [25]. ^o|| is the Ginzburg-Landau coherence length 
parallel to the sample surface at T = 0. The temperature 
dependence of ^|| is £.\\{T) = £,o\\/\/^ — T/Tc- Another su- 
perconducting parameter to be taken into account is the 
magnetic field penetration depth. A, whose temperature 
dependence is A(T) = Ao/-\/l — T/Tc, where Ao is the pen- 
etration depth at T =0. 

The Hc2±{T) curve obtained for the Nb film deposited 
on porous Si template, reported in Fig. [Ha), presents 
some peculiarities, which indicate that the superconduct- 
ing properties are influenced by the introduction of the 
porous array. In fact, if the Hc2± second derivative versus 
the temperature is plotted we can see that it changes its 
sign from positive to negative at 77 w 1.16 Tesla. This field 
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Fig. 1: Left scale: Perpendicular upper critical field Hc2± vs. 
temperature of the Nb thin film with djvt = 8.5 nm grown on 
(a) porous template and (b) non-porous reference substrate. 
The linear fits to the data close to Tc are also shown. Right 
scale: dHl^^/dT^ versus temperature. The inset shows the 
comparison between the second derivatives as functions of the 
reduced temperature of two samples, grown on the porous tem- 
plate (full circles) and on the non-porous template (open cir- 
cles), (Color online). 



value is very close to the nominal first matching field that 
we expect for the porous Si template, Hi w 1,30 Tesla, 
assuming a square porous array. This change in concav- 
ity was already reported in a previous study on the same 
kind of samples, and it was ascribed to the formation of a 
commensurate vortex structure [18], From the measured 
value of Hi it follows that the period of the porous tem- 
plate is a = 42 nm. In the following we will identify ao = 
42 nm. In Fig, [Hb) is reported the Hc2±{T) curve for the 
Nb reference film of the same thickness deposited on the 
non-porous template. As expected the Hc2±(T) behavior 
is linear over the all temperature range and the Hc2± sec- 
ond derivative versus temperature does not present any 
peculiarity except for a shallow peak near Tc. In the inset 
of Fig, [Hb), for sake of comparison, the dH^2±/d'T^ ver- 
sus the reduced temperature is reported for both the Nb 
films, in order to point out the difference in their magni- 
tude, A fit to the data close to with the expression for 



Hc2±{T) reported above, yields a value of the Ginzburg- 
Landau coherence length at T = 0, ^o|| =9,1 nm and ^o|| 
— 9.5 nm, resulting in a superconducting coherence length 
^5 — 5.8 nm and ^5 = 6.0 nm, for the Nb porous sample 
and the Nb reference film, respectively. The values of ^o|| 
are significantly smaller than the BCS coherence length of 
Nb, ^0 = 39 nm [26], indicating that our films are in dirty 
limit regime with an electron mean free path of I = 1,38 
^o|| / Co ~ 3 nm [27], Since the film dimensions in the xy 
plane are larger than ^\\{T), the expression for Hc2±{T), 
reported above, is verified in the whole temperature range. 
The Ginzburg-Landau parameter, k, — X{0) /^q|| , can be es- 
timated using the expression k = 0.72X^/1 — 9.6, where 
Al = 39 nm is the London penetration depth of Nb [26] , 
Ratios of Co||/a ~ 0.2 and A(0)/a « 2,1, measured for oq = 
42 nm, are larger than in previous works [17, 28] on perfo- 
rated Nb samples, and indicate that we are in presence of 
individual vortex pinning [29], Moreover, the pore diame- 
ter, 0, in our PS template is comparable with the vortex 
core dimension at T=0, « ^q||. This means that the 
saturation number, ns = 2|^t) ' defined as the maximum 
number of vortices that fits into a pore with diameter 0, 
is less or equal to 1, so that each pore can trap only one 
fiuxon [30], Subsequently multiquanta vortex lattice [2] 
cannot be observed in our system. 

Now we move to a more careful inspection of the R{H) 
curves of the Nb porous film. This will lead to the obser- 
vation of a peculiar behavior of these transitions, whose 
analysis represents the main subject of this work. In Figs, 
[2ja) and 2(b) R{H) curves obtained for two different val- 
ues of the temperature, T = 3,490 K and T = 3,531 K, 
respectively, are presented. 

At first glance both the curves are rather smooth and 
do not present any structures or enlargements due, for 
example, to sample inhomogeneities. However if the de- 
pendence of the first derivative dR/dH versus the applied 
magnetic field is analyzed, some distinct features can be 
observed. In particular in both the curves a small local 
maximum is present at specific values of the magnetic field. 
Let's focus on the position where the bumps, as indicated 
by an arrow in Fig. [51 start to develop. The bumps in the 
first derivative refiects the presence of a small dip in the 
corresponding magnetic field dependence of the resistance 
R{H) at the same value of H . This effect was ascribed to a 
pinning enhancement when the period of the vortex struc- 
ture is commensurate with the period of the antidots [31]. 
The bumps in the dR/dH appear indeed in our curves at 
values of the magnetic fields 7J„ when the magnetic flux 
threading each unit cell is equal to the flux quantum, $0; 
or to fractional values of $o- In Fig- H^a), where the 
R{H) measurement at T = 3.490 K is shown, the pecu- 
liarity in dR/dH is, in fact, observed at Hbump ~ 0.126 
Tesla. The period of the vortex lattice at this field value 
is a' = 128 nm, i.e. about three times the interpore spac- 
ing of this analyzed sample, ao = 42 nm. Consequently 
this field value corresponds to one-ninth of the matching 
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Fig. 2: Left scale: R{H) measurement at (a) T = 3.490 K 
and (b) T — 3.531 K. Right scale: dR/dH versus the applied 
magnetic field. In both panels the arrow indicates the field 
where the bump is present. (Color online) 



field Hi/9 w 0.129 Tesla. Similarly, in Fig. ^h) where 
the R{H) measurement at T = 3.531 K is shown, the 
bump in dR/dH develops at Hbump ~ 0.065 Tesla. The 
period of the vortex array at this field is then a" = 178 
nm, which is about four times the interpore spacing of 
this sample. Consequently this field value corresponds to 
one-sixteenth of the matching field i/i/16 « 0.072 Tesla. 
An additional bump structure is present at H ^ 0.2 Tesla. 
However, this field value does not correspond to any com- 
mensurate vortex configuration (see discussion below) and 
does not survive repeating the measurement in the same 
temperature range. Many R{H) measurements at differ- 
ent temperatures have been performed and the behavior 
of all the corresponding dR/dH curves has been analyzed. 
A selection of these curves is reported in Fig. [31 Some 
of them have been obtained by sweeping the field upward 
and downward and no hysteresis has been detected. 

For instance, the curves at T = 2.551 K and T = 3.304 
K present a bump at Hbump = Hi and Hbump — Hi /A, 
respectively. By comparison a curve with no bump, mea- 
sured at temperature T = 2.805 K, is also shown. In all 
curves the fields at which the bumps are observed are re- 
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Fig. 3: First derivatives, dR/dH, as a function of the applied 
magnetic field at different temperatures. The arrows indi- 
cate the field where the bump is present for each temperature. 
(Color online) 

lated to the first matching field through the relation: H = 
Hi/v? with 71 ~ 1,...,4. The temperatures at which bumps 
are observed, the corresponding fields and their values nor- 
malized to iJi, the values, the vortex- vortex distances, 
a, and their values normalized to ao, are summarized in 
Table [H 

We argue that the presence of the observed bumps in 
the dR/dH curves can be related to different vortex lattice 
arrangements made possible by the lattice of holes. The 
specific vortex lattice configurations occurring at the first 
matching field and at its fractional values are shown in 
Fig. 12 

In the case of Htump/ Hi = 1 a commensurate square vor- 
tex configuration is formed, where each pore is occupied 
by a fluxon and the side of this square array is just ao 
= 42 nm. Increasing the temperature the vortices diam- 
eter {k2^s) and their reciprocal distance increase, as re- 
ported in Table n When Hbump/Hi = 1/4, 1/9 and 1/16 
a square vortex lattice is again obtained with a = 87 nm, 
128 nm and 178 nm, respectively. This means that the 
pores act as an ordered template of strong pinning cen- 
ters, which is able to preserve the long range positional 
order of the flux lattice also at low fields value, i.e. at 
higher vortex spacing. As already pointed out the opti- 
mization of the vortex structures leads to the formation of 
larger square flux lattices with respect to the underlying 
artificial pinning array with the lattice constant a exactly 
equal to nao- The vortices tend to be placed as far from 
each other as possible due to the repulsive interaction be- 
tween them and at the same time they want to follow the 
imposed square potential induced by the antidots. This 
constraint gives a = aoVl^ + k^, where I and k are inte- 
ger numbers. Therefore, we should expect the fractional 
matching fields at H ^ H^/i = <f>o/a^ = *J'o/[ao(^^ + k'^)] 
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Table 1: Temperatures at which the bumps are observed, cor- 
responding fields and their values normalized to Hi, values 
at that temperature, vortex- vortex distances, a^/i, and their 
values normalized to ao ~ 42 nm. 
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0.275 
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1/16 


20.76 


178.0 


4.24 



= Hi/{P + k^) [32]. We observed bumps at fractional 
matching fields i?o/3 and Hq/4. The other bumps 

expected from the equation above at fractional fields H/^/i 
with fc ^ have not been observed. All the fields values 
at which the bumps in the dR/dH appear are shown as 
points of coordinates {H},umpiT) in Fig. (5] 
In this figure the solid lines correspond to the matching 
fields of different order, as calculated assuming an inter- 
pore spacing ao ~ 42 nm, through the formula Hi — 
$o/aQ. The dotted lines arc obtained considering a de- 
viation from the corresponding mean interpore distance 
of the order of 10% [18]. It is worth noticing that all the 
data fall into the range theoretically estimated, suggesting 
that the observed peculiarities in the R{H) curves can be 
indeed ascribed to commensurability effect between the 
porous structure of the Nb film and the vortex lattice. 
The distribution of the experimental points is consistent 
with the observation that a certain temperature depen- 
dence of the matching effect can be found for the case of 
short-range ordered templates [12]. We would also point 
out that the effect is observable in our sample only up to 
H = Hi, due to the very high value of the first matching 
field. The second matching field in fact is H2 = 2Hi — 
2.32 Tesla. From a linear extrapolation of the Hc2± curve, 
it follows that in order to see at this field a bump in the 
dR/dH we should measure a R{H) curve at T = 1.73 K, 
temperature which cannot be reached in our experimental 
setup. All the field values reported above have been calcu- 
lated assuming a square lattice. The measured field values 
do not match with the ones calculated if a triangular ar- 
ray for the pores is considered. In fact, at T = 3.490 K 
(see Fig. [D^a)) the structure in the dR/dH curve for a tri- 
angular lattice would have been observed at a field 2/^/3 
times higher than i?o/3 = 0.126 Tesla, where no peculiar 
feature has been detected. This supports our assumption 
of considering a square lattice of holes in our system. 

Conclusions. — Matching effects have been reported 
for Nb thin film grown on porous silicon. Due to the ex- 
tremely reduced values of the interpore distance the effect 
is present at fields values higher than 1 Tesla and down to 
reduced temperatures as low as t ~ 0.52. The commen- 
surability manifests both in the {H,T) phase diagram and 
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Fig. 4: Vortex lattice configurations occurring at the first 
matching field and its fractional values. Increasing the tem- 
perature the vortices diameter and their reciprocal distance 
increase, as reported in Table [1] Blue circles represent the 
holes, pink ones represent the vortices. (Color online) 
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Fig. 5; The points of coordinates {Hi,ump,T) identify the values 
of the fields and temperatures at which bumps have been ob- 
served in the dR/dH curves at fixed temperatures. The solid 
lines correspond to the different matching field orders achieved 
with the interpore spacing ao — 42 nm, while the dotted lines 
are obtained forasmuch as the regularity of the pore distance is 
achieved within the 10 percents of the average distance. (Color 
online) 
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in the R{H) transitions. The latter in particular reveal 
the formation of fractional matching states. As it was ar- 
gued in many works the vortex configuration at fractional 
matching fields are characterized by striking domain struc- 
ture and associated grain boundaries [33,34]. The pres- 
ence of multiple degenerate states with domain formation 
at the fractional field, directly observed with scanning Hall 
probe microscopy [33], seems to be high probable in our 
films. The reduced regularity of our templates, in fact, 
could be compensated by the formation of domain walls 
of different complexity. The particular domain configura- 
tion is of course a matter of energy balance between the 
cost in energy for the wall formation and the energy gain 
due to the vortex pinning. 
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